

' 'J' 


cO 

O 

00 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


<5 

o 


TECHNICAL NOTE 


oo 

rvj 

err 

o 


No. 1803 

DOWNWASH IN THE VERTICAL AND HORIZONTAL PLANES OF 
SYMMETRY BEHIND A TRIANGULAR WING IN 


SUPERSONIC FLOW 

By Harvard Lomax and Loma Sluder 

Ames Aeronautical Laboratory, 
Moffett Field, Calif. 


Washington 
January 1949 

- - 

L l_zn ^ L L L 

_ i- ", 

L ^Ll ur: An j L 




f a ; 



UJ 

D3 -~- — 



V , J 



TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 


HATIOHAL ADVISOR! COMMITTEE TOR AEROMUTIU, 


lllllllll 


DIMM' 


TECHNICAL SOTE HO. 1803 
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SYMMETRY BEHUTO A TRIAHGULAR WING HT 
SUEERSOHIC FLOW 

By Harvard Lomax and Loma Sluder 


SUMMARY 

A method developed in a previous report for finding the induced 
velocity field behind a supersonio wing with known, load distribution, 
is used to find the downwash behind a triangular ving with subsonic 
leading edges. Results are given for the chord plane in the extended 
vortex wake of the wing and for the vertical plane of symmetry up to 
about 20 percent of the wing span above the plane of the wing. 


IHTROEfUCTIOH 

A method has been developed in reference 1 for computing the 
downwash behind wings of known loading flying at supersonic BpeedB. 
The solution was based on the distribution of supersonic doublets 
over the plan form and wake of the wing in a manner determined by 
the load distribution. The method was applied in reference 1 to 
the calculation of the downwash behind a triangular wing with 
leading edges swept behind the Mach cone from the vertex. Hear the 
trailing edge, however, these calculations were limited to a region 
close to an extension of the center line of the wing, that is, the 
x axis (fig. l), and in. fact were exact only on this line itself. 

A simple first approximation was advanced for the downwash variation 
about the x axis j namely, that the difference in the downwash at the 
position (x,y,z) from the value at the position (x,0,0) was a linear 
function of the z distance and independent of the y distance. 
However, the accuracy of this approximation could be tested only for 
large values of x where the downwash was known for all values of , 
y and z. 

The purpose of t his report is threefold: First, to continue 

the exact calculations so as to include all points on the xy and 
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xz planes behl nd the trailing edge for distances up to about a 
semispan from the x axis in the xy plane and up to about 40 percent 
of a semispan in the xz plane; second, to compare these values of 
downwash with those obtained from the approximate method given in 
reference 1; and, third, to serve as a guide through some of the 
more difficult mathematical manipulations so that the calculations 
can be extended to other plan forms. 
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LIST OF IMPORTANT SIMBOIfi 
velocity of sound in the free stream 
root chord of w ing 

couplets elliptic integral of the second kinfl with 

modulus k,}^, respectively 

incomplete elliptic integral of the second klnfl with 

«t 


argument t and modulus k E(t,k) = J ' dt 
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complete elliptic integral, of the first kind with 
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incomplete elliptic integral of the first kind with 
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V(2^i) 2 -P 2 (y-yi) 2 -0 2 (z-2:i) 2 


VL,V,V 


perturbation velocity components in tbe direction of the 
x, j, and z axes, respectively 


AUg Uu-Uj 

Y c free-stream velocity 


YP 


z component of velocity induced by doublet, distribution 
over plan form 




z component of velocity induced by doublet distribution 
over wake 
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Cartesian coordinates of an arbitrary point 
Cartesian coordinates of source or doublet position 
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angle of attack 

V«o 2 -l 
P tan 

density in free stream 
perturbation velocity potential 

semi vert ex angle of triangular wing 
sign denoting finite part of integral 
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Subscripts 

u 

l 


L.E. 


T.E. 


¥ 

P 

s 


A,B,C,D,E 


conditions on 
conditions on 
conditions at 
conditions at 
wake 

plan form 

conditions on 

conditions in 
wing (figs. 1 


upper portion of surface 
lower portion of surface 
leading edge 
trailing edge 


discontinuity surface (at z^-O) 

regions A,B,C/D, or E, in wake of triangular 
and 2) 


THEORY 

The theoretical development in this report is subdivided into 
two parts. First, the problem of determining the downwash behind a 
flat— plate triangular wing swept back of the Mach cone is reduced to 
the solution of certain types of elliptic integrals. This material 
has already been presented in reference 1, and. the brief description 
presented here should be a sufficient review. Second, an evaluation 
of the resulting integrals is made in both the xy and xz planes. 
Some of the detailed evaluations are given because the methods and 
substitutions presented may be useful in the solution of downwash 
problems for other plan forms. 

Boundary Conditions and Their Solution 

The load distribution for a triangular flat-plate wing (fig. l) 
has been determined und er the assumptions of t hin— airfoil theory and 
Is given as 


An _ 2Aus _ ^ 0 2 axi 

* V ° 3,3 VeoW-P^x 2 


( 1 ) 


Since the loading coefficient determines the jump in the u induced 
velocity in the plane of the wing, the definition of the perturbation 
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velocity potential 4> gives 


A<t s 



Aug dx 


where ££> B represents the jump in the velocity potential in the xy 
plane. In any lifting-surface problem there are three regions in the 
xy plane for which the equations for A$ s are different: the region 
defined by the boundaries of the lifting surface itself, the region of 
the vortex wake behind the lifting surface, and the remainder of the 
plane. For the triangular wing the jump in potential over the lifting 
surface is given by the equation 

A$ s = H V0 o 2 2i 2 -0 2 7i 2 , (2) 

where 

x i £ c o 

and 

E =. 22a 

E 0 e 


In the vortex wake (the semi— infinite strip behind the wing extending 
from tip to tip with sides parallel to the free— stream direction), 

A<P S is independent of x x and is given by 

. A% = H /\/e a 2 c 0 2 -p 2 y 1 2 (3) 

/ 

where 

xi > c Q 


Everywhere else in the xy plane A2> s is zero. 

The basic partial differential equation satisfied by the pertur- 
bation velocity potential in supersonic flow is well known to be 


dx 2 dy dz 2 


w 


The general solution of equation (4) is given in reference 2 in the 
form 
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-m b 



dx* dya. 


( 5 ) 


where 

r c = A /(z-^i) 2 -P 2 (y-yi) 2 -P 2 (z-z 1 ) 2 

and the subscript s indicates that the function is to be evaluated 
at Zi equals zero. The region t is the portion of the combined 
area of the plan form and wake "bounded by the leading edge of the 
wing, the sides of the vortex wake and the trace in the Zi=0 pl ane 
of the Mach forecone with vertex at the point (x,y,z). The Bign f 
is to be read "finite part of" and is used in both references 1 and 2. 
It has the property that 




f(x)-f(b) _ 2f(b) 

0 «> 3/ * 


(6) 


Since the particular problem of this report is a lifting plate 
without thickness 


d4»u _ 

3zi 8zi 

and equation (5) reduces to the form 
t(x,y,z)= t(x,y, z) p + t>(x,y,z) w 




plan form 


AQ a dxi dyi 

[(x-x 3 .) 2 -p 2 (y t -yi) 2 -f 2 z 2 ] ^ 


z|3 2 

2rr 


f [' A 1 > h dx x dyi 

•- [(x-xi) 2 -P 2 (y-yi) -P 2 z 2 ] 


( 7 ) 


Solution in the xy Plane 

lr i g convenient to divide the study of downvash in+o t : parts; 
one , the study of the offers of the doublets distributed ov ~ <he 
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plan form; and, two, the study of the effects of the doublets distri- 
buted over the wake. Such a division has already been implied by the 
arrangement of equation (7). Hence, In the following and Vp 

will be derived separately and the total dcwnwash will be their sum, 
w = Wfr + Wp. 

Effect of doublets in the plan form .— For the purpose of integration, 
it is convenient to divide the area behind the wing into three regions 
as shown in figure 1. The division line separating these regions are 
formed by the Mach cone traces from the traj l in g-edge tips. The limits 
of integration which form the bases of the division differ in each of 
the three regions and are discussed in more detail in reference 1. 

The following symbols will, be used in the derivation of the 
expressions for downwash, in the xy plane, induced by the distri- 
bution of doublets over the plan form of a triangular wing swept 
behind the Mach cone: 

Ex, E 2 , Eg complete elliptic integrals of the second kind with 
moduli- ki, k 2 , and , respectively 

Ki, K 2 , Kg complete elliptic integrals of the first kind with 
moduli ki, k 2 , and kg, respectively 



g(p-fp , )-V2g(p t -p) (p+g) (pMQ 

(jL-p , +25 

(pp»+i 2 )+ V(i g -p 2 ) U 2 -^ t2 ) 

p+p i 
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T1 




(^ , +g 2 )+ A /(g 2 -ix 2 ) (S 2 nx ,a ) 

§(^+}i»)+V2g(^ t -^) (n+5) (p'-T) 

(ti|i t +i 2 )-V(i 2 -ii 2 ) ( S 2 — i-l t ^T 


P(y-yi) 

Co ' 



P* 


*o 


+ 2fl*i 

Co 


I 

Co 


to Xq -1 


The downwash wp(x:,y,0) may "be obtained "by considering lim 

jk> 0 dz 

in equation (7). It can be shown that, in this case, this limiting 
process corresponds to taking the partial derivative of equation (j) 
with respect to z and then simply setting z equal to zero. Thus 
the expressions for Wp in the regions A, B, and C are, respectively. 


HP 

wp„ — 

r A 2ir 


A 

J to W » 


'■L 


^ VCu’-J]) (if-p) 


(ivr 


dri 


= -M 
2 * 
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The solution, of the three integrals Ii, Ip, and Iq, will he lie cussed in Appendix. A. 

The expressions for downvash in regions A, B, and C Bay then he expressed as the following 
b ingle Integrals which can he handled by standard numerical methods: 
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Effect of Doublets In the Wake .— The study of the downwash induced 
on the zy ' plane by the doublets distributed over the wake will also 
be divided into the three regions indicated in figure 1, and the symbols 
listed as follows will be used in the derivations: 


a A 


a B 


ac 


&A 

?A 



7b V&b V V 
7c (v'- & c ') 

S C W v ‘^ 


E C 



conplete elliptic integrals of the second kind with 
moduli k^, kg, and kg, respectively 



incomplete elliptic integrals of the second kind with . 
arguments l/&k, l/ap, and ac/kc, and with moduli k^, 
ks, and kc, respectively 


k A’ E B’ *C 



couplets elliptic integrals of the first kind with 
moduli k&, kjj, and kc, respectively 



ineamplete elliptic integrals of the first kind with 
arguments l/a^, l/ap, and a^/kg, and moduli k^, kg, 
and, kg,- respectively 
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(w»+i 0 a )— V , (i 0 2 -v 2 ) (S 0 2 -v» 2 ) 
v+v f 

g 0 (v+vM-'/2g 0 (v^v 8 ) (v-6q) (lo+vQ 
v'-v+2£ 0 

(w , -M 0 2 )WU 0 2 -v 2 ) (Sq 2 -*' 2 ) 

v+v* 

(vv , +g 0 2 )+V(So 2 -^ 2 ) (io 2 -v ?2 ) 
v+v ? 

i o (v+V«)H-//2g o (v-V t ) (V-l 0 ) ( Ip+V'M 
V*-V+2£ q 

(vv*-0-^ 0 2 -v 2 ) (lo 2 -v^ 
v+v* 
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n 

P(y-7i 

c o 

V 

7 o+ 8 o 

V* 

Ho - ®o 

£0 

z o“ 1 


On integrating the second term in equation ( 7 ) with respect to 
and using the notation just presented 


0 — 5E ^o z o c o 

¥ 


(n-v 1 ) 


2lt P J 1-2 ( T l 2+Z 0 2 ) V to 2 — 1 \ s — z o 2 


dt) 


(HO 


The limits on the integral as previously noted differ in the 
regions A, B, and C shown in figure 1; however , in each case the 
limits are roots of one of the two radicals in the integrand. 

It is desirable to express equation (l4) in a different form in 
order to obtain an expression for downwash in the plane of the 
airfoil. Integrating by parts 


HP c Q 
2it p 
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J (v— T]) (tj — v ’) tan 
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z OrJ I O 2- T l 2- z o 2 




V+V>-2T1 

2/s/(v-n) Ol-vO 


tan' 


— 1 


So 1 * 


z oJio -T[ Z -Zo Z 
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When 11m — Is considered, It can "be shown that In all three regions the contribution to $*• 

z 0 -M) ok 

the dowrurash made by the doublets in the wake Is given by the expression 


% „JL. r Ll - V+Vt -^1 JT^ an 

^ Jjc a lV(^)(Tl-Vf) * ° 


(15) 


Ihe solution of equation (l?) in regions A, B, and C will be considered separately. 

Region A 

In region A, Li«v and Ls«v*. die substitution t) » ^ eliminates the i -f nsn-r term in 


1 +t 


the radical of the integrand, and equation ( 15 ) becomes 
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5 A“ vT 


v+v* 

r 11 
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Integrating the second term "by parts and applying the fundamental properties of even nnfl 
odd functions yield 



ntiR Jacobian transformation cu = sn u reduces the integrals In equation (17) to standard 
elliptic forms (reference 3), and 
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Region B 


In region B Lx = £ and L2 =■ v ’ and equation (15 ) may "be written 
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*rrn f ^ 
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2 *5, 
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a/[ (V— TO a o +T))l [Uq-TT) (H-V 1 )] 


TB+Sgt &g— V T 

The transformations ^ = — and (a = r t where 
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and 


(£ q + 7b) ( v_6 b) + +e b) (v-7b^ “ 0 


reduce equation (19) to 
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(v+V l ) 7 B(l +a R ai ) 
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VVX-^-Tar-; 
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and, after algebralo simplification, may Be written 
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In region C f Li = 5 and L2 » — | . and equation (15) is written 


v w, - 


S3 


C 2 utg 


-/°- 

2i 


V+V*— 2q 


2Ti*/(v-n) (n-v) 






( 22 ) 


The 


derivation of the expression for w^ c is similar to that for vy A with the exception 
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Solution in the xz Plane 


Just as in the study of downwash in the xy plane, so also in its study in the xz plane it 
is convenient to consider separately the effects of the doublets distributed over the plan form 
nnfl the vahe. The subscript notation for w^r and vp is the same as before and again v is equal 
to v w + Vp. 

Effect of doublets in the plan form .— In the xz plane two regions are indicated in figure 2. 
Eeglon E lies between the Mach vedge from the trailing edge and the line of intersection of the 
two cones from the trailing— edge tips. Eeglon p connects this region to infinity. Again the md 
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limits of integration form the basis of the division into the two 
regions . 

The symbols listed below will he u^ed in the derivation of the 
expressions for downwash in the xz plane , Induced by the distribution 
of doublets over the pi an form of the wing. 


complete elliptic integrals of the second kind 


E 4' E S' E 4 * E 5 . , 

0 0 with moduli k 4J kg, , and k^, respectively 


I. XX complete elliptic integrals of the first kind 

0 0 with moduli k 4 , k 5J k 4 3 and k 5 , respectively 

0 o 


k 4 


fj (x— Xi ) 2 — P 2 Z 2 


k 4o 




^(x-Xx ) 2 -^ 2 

So 1 ! 

VC^Gp )^^ 5 

6q°o 



Region D 


In region D, the first, term In equation ( 7 ) Is written 
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1 ^ y°o -p yi 

/ / 7: ,£ & a k 27072 ayi 

l(2MCi) -P yi -0 a J 


(24) 


Integrating with rasnect to vi 




0 o x l 


(K^-E^) dxx 


C h anging the variable of Integration 


. _ *HP T ’° 0 o L xi 4 -Ve n 2 x 2 +P 2 i £ (k 4 2 -e o 2 )_ / K4-^ + \ a . 


Taki n g the partial derivative of 4 p with respeot to z gives the expression for the downwosh as 
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0 o 2 °o 2 V^ 4 O 2 - 0 O 2 / V^oWa 1 2 (k 4 0 2 -^o^) 

/ Ai 0 o g -A / go \ 

\[0 o 2 x 2 +p a z 2 (k +o 2 -0 o 2 ) ] a / 2 / \ 3c4 2 -g 0 2 / 


|*4 0 - V 0 O V+^z 2 (ti o 2 -0 o 2 ) J 

(K 4 ^ 4 ) dk 4 j. 


Region E 

In region E, a similar derivation yields the solution 
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C\g e 0 c 


k 5 l-k 5 e 0 o e Ui 2 k 5 2 ©o 2 +^ 2 ^ 2 (l~k5 2 0o 2 ) ] S/,S 


-i} 


’ f (K4-E4) 


x^ffo 2 


g 0 r 

k4 2 -flo 2 l 1 0o 2 x 2 +p 2 z 2 (k 4 2 -^o 2 ) ] 3 / 2 




dk* 


(26) 


Effect of double ts, In the Treks .— The limits of integration again necessitate the division, of the 
portion of the plane behind the trailing— edge wave into two regions D and E. Tdw following list of 
symbols will be used in the derivation of dnwnwash in the xs plane, induced by the doublets distributed 
over the vortex wake: 
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a D 

V(x-c 0 ) -p z 
Pz 


a *D 

A /(x-c 0 ) 2 -p 2 z 2 

(x-c 0 ) 


ag 

e o c o 

pz 


a ’E 

^o c o 

7p 2 z 2 ^o 2 Co 2 


Ed, Eg 

complete elliptic integrals of the second kind vith 
moduli kj) and kE, respectively 

E ( a V k 'D) 

E(a»E,k'E) 

j 

incomplete elliptic integrals of the second kind vith 
arguments a' E and a' E and moduli k' E and k* E; 
respectively 

Kd> Ke 

complete elliptic integrals of the first kind vith 
moduli kpj and kE, respectively 

T 

E(a«D,k’D) 

F(a* E ,k' E ) 

incomplete elliptic integrals of the first kind vith 
f arguments a’p and a*E and moduli k'pj and k*E. 

respectively 

k D 

®o c o 

c 0 ) 2 -p 2 z 2 


k *D 

Vi-^d 2 


k E 

y(x-c 0 ) 2 -P 2 z 2 

0 o c o 


k *E 

' yi-^E 2 




sn* 


! (ic sq) 


“** a D 2 


B^ 2 f *f n f — 'l 


“ifi 


Region D 

In region D, the second term in equation (7) becomes 


*w = - 

*D 


»Tro2 


^o c o 

r p 


i-PV7i E+ zS 
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/ Ve ‘c 0 ‘-^>i‘ d 7l J 
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I^Ti 


[ ( X-Xi ) 2 -p 2 ( 7i a +a S ) ] 8/E 


Integrating with respeot to x x , and. using the definition of 

*JQo^ c o Z -^Ti Z 




d o°o 
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(y x e +a 2 ) V (x-c 0 ) E -P 2 yi 2 -P 2 * 2 


dy x 


H(x-c q ) ^0 o 0 r Jl-^aF- 

fl d_ . // _ _ \e q'e^e J 


PzV(x-o 0 ) 2 -^?^ f 


2^2 
•D 03 
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ro 
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where 
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Applying the Jacobian transformation - sn u 

4 > - E(x-c 0 )e 0 z o 0 s l~En Z u du (28) 

W D *fi Z J(x-c 0 ) s -fi 2 z *4> 

The expression for 4u as given in equation (28), is integrable (reference 3 ) and becomes 

W D 



5(x-c 0 )6 0 2 c q 2 
*{3z J{x-o 0 f~^z^ 






- i^a^k^) - U(a» D ,k» D ) 


5d11 

k d JJ 




t? k d 
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a/i^'d^'d 2 

*Jl-a 'p" ' 


% P(a» D ,k» D )^ D F(a l D ,k' D )-K D E(a» D ,k* 



(29) 
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Since it may "be Bhown that 
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[Eq F(a» D ,k* D ) - ISjj F(a» D ,k T D ) - K D E(a» D ,k» D ) ] 


> e ^ 0 % 2 L 


(KjrSp) a'xrKi) 


(l-aVW 


the ejcpreaaion for the doynyaah obtained "by taking the partial derivative of with respect 

to z may "be written 


HP r to. + K s( :L - a 'P a ) , _ ^5pL [K B r(a> D ,k«D) - % P(a' B ,k* D ) - K D B(a* B ,k» D )]} (30) 

WD * 1 a'u a*D J 

Region E 

In. region E a Bimilar derivation gives 

(“I? ~ a 's 2 %) - Vl-a’/ 't% Pta'j.k'j) - Ej P(a> E ,k> E ) 

(31) 
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DISCUSSION 

In reference 1, the downwash behind the triangular wing was 
computed along the x axis only; hut it was indicated that only 
slight variations with y should he expected for some distance on 
either side of the axis. The value of w/Wq 1 along the x axis 
derived in reference 1 for various 0 o ’s is shown in figure 3* 

The results presented in figure 4 can he used to assess the 
accuracy of the prediction of reference 1 in regards to the spanr- 
wise variation of downwash. The downwash in the xy plane is pre- 
sented for various 0 o ’s and spanwise stations and for all values 
of x from the trailing edge to a point where the asymptotic value 
is closely approached. In figure 4 the curve for y Q = 0 represents 
the results of reference 1. The region covered in the y direction 
extends from the x axis out to about i 0„, where in the coordinate 

2 

system used x Q equals x/c Q , y Q equals 3y/c Q , a n d 9 Q is the 
semispan of the wing. Within about a half span from the trailing edge of 
of the w ing no general statement can he made as to the variation of 
w/w Q in the y direction except that in the region considered the 
increase or decrease in the values of downwash from those at y Q = 0 
was less than 0.1. For distances greater than a half span from the 
trailing edge, however, the variation is quite uniform and w/w Q 

deviates from its value at y=0 only slightly for - ^ e 0 <J 0 < \ e o‘ 

0 

For 0 o =O.6, a more extensive Btudy was made of the variation 
of downwash with y. Figure 5 presents values of w/w D across the 
span for several positions behind the trailing edge. Immediately 
behind the trailing edge the value of w/w Q falls and approaches 
-00 as the wing tip is approached. However, at 0.4 of a root chord 
behind the trailing edge (x Q =1.4), w/w Q riBes and reaches the 

value of 0.7 as the wing tip is reached. At x Q =2.2 the spanwise 
variation of w/w c is essentially constant. The approximation of 
reference 1, that the downwash is independent of y in the neigh- 
borhood of the x axis. Is also shown in figure 5- It i 0 evident 
that this approximation is useful out to about a third of a semispan. 


1 The variable w/w Q (i.e., (wp + w w ) /w Q ) represents the total 

downwash behind the wing divided by the induced vertical velocity 
on the wing itself. If e is the downwash angle and a the 
angle of attach of the wing then 

w _ de 
w c “ da 
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The variation of downwash in the xz plane is presented in 
figure 6. The curves represent values of v/w Q from the trailing— 
edge wave downstream to a point where the as ympt otic value is closely- 
approached. In the immediate vicinity of the Mach cones from the 
trail in g-edge tips (i.e., x Q SJ 1 + 0 O ) the curves were not continued 
because w/w Q becomes very large and approaches negative infinity 
as the Mach cone is approached. (Since this effect results from 
infinitely large values of the radial component of induced velocity 
at the Mach cone, it does not exist in the z Q =0 plane.) Such a 
behavior is consistent with .the mathematical idealization of infinite 
pressures at the leading edge and of an abrupt fall of load at the 
trailing edge. However, in an actual flow field where these 
phenomena do not exist the flow will experience a milder change in 
passing across the Mach cone. Even in the theoretical results 
presented in this report the growth of the vertical induced velocity 
in the neighborhood of the Mach cone is logarithmic, and the interval 
in which w/w Q is appreciably distorted from the general trend is 
very small. 

Some further insight into the behavior of w in the vicinity 
of the Mach cone from the trail ing-edge tips can be obtained by 
studying a single vortex which extends infinitely far ahead from the 
origin at an oblique angle to the flow and infinitely far behind the 
origin parallel to the flow (fig. 7)- The half of the vortex which 
extends ahead makes an angle with the free— stream direction lesB than 
the Mach angle so that the component of free-etream velocity normal 
to it will be subsonic. Thus, outside of the Mach cone originating 
at the sudden bend in the vortex at the origin, the flow will be 
exactly likB that of a linearized compressible subsonic vortex with 
a super Imposed uniform velocity parallel to the line of the vortex. 
Inside the Mach cone, however, the flow is completely changed. 

Figure J gives an indication of the change. The term "bent" vortex 
refers to the vortex along the x axis which is turned suddenly at 
the origin from the angle it had maintained from —m. The term 
"unbent" vortex on the other hand refers to a vortex which maintains 
the same angle from — » to +“ . Tbe unbent vortex is included 
to show the effect of the sudden turn. The figure shows that on 
the z=0 plane (section AA) the downwash is finite and continuous 
in passing through the Mach cone, but that above the z=0 plane' 
(section EB) the value of w becomes infinite as the cone surface 
is approached from the inside. This behavior at the Mach cone 
may aid In interpreting the discontinuity in the results for the 
conplete wing as given in figure 6. 

A conparison of the simplified study made in reference 1 with 
the results of this report in the xz plane is given In figure 8. 
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The conclusion of reference 1 was that to a first order the variation 
of w/w Q with z o /0 o was linear with a slope -(l/E 0 ).. The figure 
shows that the approximation is useful up to ahout a third of a semi- 
span. 

Values of v/w 0 were not computed for points off of the xz 
and xy planes although the methods employed are quite capable of 
handling the problem. The results already given would apparently 
indi cate that the approximations of reference 1 are valid in the 
vicinity of ~0-/3)&o about the x axis. This assumption can easily 
be checked for large distances behind the trailing edge by consider- 
ing the flow field as Xq approaches Thus figure 9 shows a com- 

parison between the exact value of w/wo derived by means of the 
linearized equation and the approximate value given in reference 1 
for all points greater than about one chord length behind the 
trai ling edge, Once again the agreement is fairly good out to about 
one-third of a semispan either vertically or horizontally from the 
x axis. 

In order that some idea of the variation of downwash with Mach 
number at various positions downstream could be obtained, "figure 10 
was prepared. This figure shows values of w/w D on the x axis 
plotted as a function of Mo for various values of Xq . The value 
of the sweepback angle is 45° and the Mach number range covered 
could be extended to 1.4 and the leading edge of the wing would still 
be subsonic. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 


AFEENDIX, A 

EVALUATION OF SEECIAL INTEGRALS 


Integral I x 


Since there acre no singularities in I x , the finite part sign 
may be discarded. The linear term in the radical is eliminated by 
the transformation t] = ( 7 x +S x t)/(l+t) and the integral becomes 



< jO 

O 


'(u'-7i) (7ini) 

■5" 


~ 8 1“U 


(8i^7i)7i H , -7i 


£ 


2 

t 2 

1 - 


3/2 


dt 


(Al) 




The expressions for 6* and 7i may "be combined to give the useful identities 

% = 7i&i 


(7i-u) = (tiM'i) 


Hotlng that the Integrand is an even function, equation (Al) may be reduced to the 
canonical form 


by the substitution 



Vl -< o g 

(l-fci E <n s ) 8 / 2 


dm 


(A2) 


By the introduction of the Jacobian elliptic functions (reference 3) in the transformat ion 
tD =» sn u, the integration may be completed, and 
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■ ■■ ■ 


(A3) 


where od u s on u/dn u. 


Integral I 2 


Aq the first step in reducing 1^ to canonical foiro the integral is written 


la 



■. , ! 2 “£(h+h' )+wi’ 

£ 2 +£(h+h' )+ pn« 


L 25(«) 

2£(i H ) -l 


Jtl 


VU+n) VU-n) (thO 


In this case the following identities may he obtained directly from the definitions of 
7z and B 2 . 


m 


Cr JJ -u)(^-B 2 )+(r 2 -^)(n-8 0 ) = 0 


and 


(i+7 2 )(^-6 2 )+((i*-7 2 )(|+5 2 ) o 0 
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By applying the fundamental properties of even and, odd. functions, the first two integrals 
in equation ( A5J can readily he integrated. The procedure for handling the finite part Blgn 
over the third integral xill he considered in detail. Since, 




NACA TB Bo. 1803 



then "by equation (6) 


1 

r dm n 

_ 1 1 q 

p A m r dflD 1 

^0 ( l-KD 2 ) 0/2 

To (l-ci 2 ) 3/ ^Vl-^ 2 J o 2° /C (1 ^d)' J; '/l-hs 2 72 ( 1 — k e 2 ) J 


= 2 


I 


I j/J) n | . | 

i _ Vsu-®) (!-**")•{ 


:i-ife 2 )] 


J1 2 


r du _ 1 . i 
_J OT^u ^l /V /2(l-oi) (l-k£ 2 )- 


f f nj a/l-k^to 2 S2 1 _ 1 . 1 

2 il*® + (i-kgS) “ l-k/J CD .’ >1 Vs(l-a>) (l-ka s V 




The solution of equation (A5) becomes, after algebraic simplification, 

' a a (n-ii»+ 2 £)-gn£l - — [ZESJES. S2 l 

t £ J 2 £ 2 V 5 l 2 J 


-{/ 


k2 


(Hi) (Hi 1 ) 


(A6) 


LO 

l_U 
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Integral Ig 


The procedure for integrating I 3 is similar to that previously 
discussed in connection with lx. In this case, the integral is 


canonicalized hy the substitution 
written 



and the solution may be 


Is » 



(A7) 
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Figure 3. - Variation of the y downwash on the x axis downstream from the trailing 
edge (Reference t ) . 
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Figure 5. - Variation of downwash across the span at various 
stations downstream > of wing trailing edge for 0 O = 0.6. 





Figure 6. - Variation of 
from the trailing edge. 













figure 7 - Induced vertical velocity field for bent and unbent supersonic oblique vortex 
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(a) 6 0 -- 0.4. (b) & 0 -- O. 6. (c) 0 O = 0.8. 


Figure 8. - Variation of downwash in the x Q z 0 pfane at various positions on x 0 axis. 
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Figure 9. - Downwash at a targe distance 
behind triangular wing. 





